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The paper discusses mechanisms for decay of supercurrents in ultrathin superconducting wires 
driven by quantum fluctuations. We argue that momentum conservation strongly suppresses proba- 
bility of such decay and estimate the rates for two decay channels: potential scattering of condensate 
due to the disorder and attenuation of the plasmon mode due to the presence of normal component, 
i.e., Ohmic losses. We find that while both mechanisms yield non-zero decay rates, their values are 
too small to provide any substantial contribution to the resistivity of the wires. The rate associated 
with the latter mechanism, however, is much greater, and it is possible that under the appropriate 
conditions dissipation may lead to appreciable enhancement of quantum phase slip transitions. 
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Properties of superconducting nanowires at ultralow 
temperatures have been a subject of extensive theoreti- 
cal and experimental studies for several decades. It is be- 
lieved that such properties are controlled by proliferation 
of the so-called quantum phase slips (QPS), - topological 
fluctuations resulting in discontinuous changes of phase 
of the superconducting order parameter driven by quan- 
tum fluctuations [ij . Such QPS are expected to lead to fi- 
nite resistivity of superconducting wires at temperatures 
much lower than critical temperature (Tc) and, under 
certain conditions, to superconductor- insulator quantum 
phase transition [l|-l3|. 

Despite theoretical predictions that QPS should be 
observed in sufficiently thin nanowires [2|, their exper- 
imental observation turned out to be rather difficult and 
inconclusive. While initial experiments seem to have 
demonstrated the existence of QPS fluctuations in MoGe 
nanowires subsequent measurements did not con- 
firm these observations [Q]. More recent measurements, 
though, have suggested that QPS might have been ob- 
served in nanowires carrying sufficiently high bias cur- 
rent, close to the value of critical current (6;J. 

In this paper we estimate the rate of the QPS in su- 
perconducting nanowires. We find that this rate is con- 
trolled by two factors: (1) semiclassical exponent due to 
the quantum tunneling of the order parameter through an 
effective energy barrier, which arises due to the suppres- 
sion of the superconducting density at the QPS cores; (2) 
a probability for the superflow to change its momentum 
by hp/2, where p is 1-dimensional (ID) density of con- 
duction electrons. We argue that in experimentally ac- 
cessible regimes the latter factor is responsible for major 
suppression of the QPS rates. Particularly we estimate 
the rate for two mechanisms of momentum relaxation, 
such as scatteringof the condensate at the disorder po- 
tential, e.g. Ref. Q , and due to the dissipation of the 
plasmon mode resulting from its coupling to the normal 
component. We find that potential scattering mechanism 
is ineffective in wires with sufficiently large number of 
transverse channels (i.e., whose radii significantly exceed 



the interelectron distance), which corresponds to practi- 
cally all up-to-date QPS experiments. The dissipation 
based mechanism leads to much higher QPS probability 
and, may possibly be an explanation of experimental ob- 
servation of QPS events at sufficiently high values of bias 
currents [H]. 

Dynamics of low energy excitations of a superconduct- 
ing wire is described by the following Lagrangian density 
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[(a.0)2 + c2(9.0)2], (1) 



where (/)(x, t) is the phase of the superconducting order 
parameter. The variables x and r are the coordinate 
along the wire and the imaginary (Matsubara) time re- 
spectively. Effective Lagrangian in Eq. JT]) can be de- 
rived from microscopies, e.g. Refs. 0,0]) as an expan- 
sion in powers of gradients of (f>. The last two terms 
describe propagating plasma mode [1]; here C is capaci- 
tance (per unit length) of the wire and Cs is the (Mooij- 
Schon) phase velocity Q. Note that the plasma mode 
in Eq.(IT|) is gapless, which is a property of ID wires: In 
sufficiently thin wires the screening turns out to be effec- 
tively weak because the electric field is "pushed out" to 
the region outside the wire. Note that Eq.dT]) is valid only 
at sufficiently large distances (and times), greater than 
some cut-off length ^, of the order of the superconduct- 
ing coherence length. The term ihpdr4>/'2 is the so-called 
Berry phase [lO|. The iniportance of this term has been 
discussed in Ref. 0, 0, B Ei : It accounts for momen- 
tum conservation during the phase slip formation, which 
leads to the suppression of the QPS rates (see below). 
The term proportional to dx4> is due to the applied bias 
current /. 

We start by considering a translationally invariant sys- 
tem. The rate of the QPS events can be evaluated by 
using the instanton method We are interested in 

the decay probability of a current-carrying state |/), i.e., 
with bias supercurrent /. A corresponding instanton so- 
lution, i.e., a classical trajectory for the Lagrangian in 
Eq. ([T]) satisfying appropriate boundary conditions 
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consists of a kink-antikink pair, 

(j)c{x, t) = tan"^ — ^ — - tan"^ — , (2) 
and the rate is given by the expression 

rqps =Klm J d{Ti - T2) J dxidx2e~^<'^^^^'' . (3) 

The prefactor K in Eq. ^ is related to the so-cahed 
fluctuation determinant 12], whose value can only be 
estimated, e.g. Ref. K ^ c^/^^. The classical action, 
corresponding to the classical (instanton) trajectory from 
Eq. (12), is readily obtained as S{(f>c) = J drdx Co{(j)c), 
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In Eq. (gl a = nhCcs/'ie'^ and Aa; 
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Ti — T2 . It is convenient to carry out the Ar integration 
in Eq. ([2]) and then do the remaining Ax integral. The 
Ar integral is formally divergent and one first needs to 
analytically continue it by bending the contour of inte- 
gration into the complex plane in the region of negative 
argument 1^ . As a result one obtains 
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The subsequent integral over Ax yields 0. To see this 
one may use identity (l3| 
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Thus, replacing the Bessel function in Eqs. ([5]) by its 
integral representation, Eq. ([6]), the integral over Aa; in 
Eq. ([3|) yields (5(ecsp± Jcos6'), and therefore the remain- 
ing integral over 9 is nonzero only for / > ecgp. It is easy 
to verify that for a superconducting wire quantity ecsp 
always greatly exceeds the value of critical current [6| and 
therefore Fqps — 0. That is, the momentum released as 
a result of a phase slip event needs to be absorbed by the 
plasmons. Such process, however, is suppressed by the 
Landau criterion i.e., impossibility to transfer both 



momentum and energy from a superfiow to a bath of ex- 
citations (i.e., the plasmons), whose spectrum is linear, 
e.g. Eq.©. 

If translational invariance of the system is broken, 
the momentum conservation prohibiting QPS formation 
is "violated" and, as a result, the QPS rate becomes 
nonzero [ll[ . The presence of lattice or disorder (in 
the following we will study the effects of the latter) ob- 
viously modifies the Lagrangian density in Eq. ([T]). In 
the second order terms the disorder leads to the renor- 
malization of the stiffness of the plasmon mode, e.g. Ref. 



Q. Such renormalization, i.e., averaging over the dis- 
order realizations, is well justified since the disorder po- 
tential varies on a scale small compared to cut-off length 
^. The contribution from the Berry phase term, how- 
ever, comes from the discontinuity of the phase at the 
discrete set of points, e.g., QPS centers, as we have seen 
above. Therefore the electron density p in this term does 
not self-average, i.e., one should explicitly account for its 
spacial variation. Since the classical equation of motion 
is independent of the Berry phase term, the instanton 
solution and the rate are still given by Eqs. (I21E1), while 
the classical action becomes 
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The only distinction between Eq. ([7]) and Eq. (U) is 
obviously the last term. By writing p{x) = po + ^p{x) 
and assuming that po ^ 5p (which is a well justified 
assumption for a wire whose diameter is much greater 
than the interelectronic distance, - recall that p is ID 
density, i.e., averaged over the wire's cross-section), we 
rewrite the x integration in Eq. ([3]) as 



dx^dx^e''' "'"'^''^Fix, - X2) ^ {L/pD 
X / dAx{Sp{Ax)6p{0))F{Ax)e'''P''^'' , 
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where F is the remaining x-dependent part of the inte- 
grand function in Eq. ([3]). In deriving Eq. (jSl) (together 
with the inequality po ^ Sp) we used a fact that F{Ax) 
varies in a scale Cs/I, which is much greater than p^^. 

Then, substituting Eqs. ([51 [H]) into Eq. (jS]) and using 
Eq. ([6]), after some straightforward algebra, one obtains 
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where (|(5ppp) = J dx{6p{x)6p{0))e^P^ and po = t^pq. 

The first two factors in the rhs of Eq. ^ can be put 
in a traditional form ujQe~^°, where ujq ^ CgL/^"^ is an 
attempt frequency and = (2a — 1) In (cCs/ttIS,)- This 
result has essentially been obtained in Ref. [2| . The third 
factor (which we will denote P3 below) is a correction due 
to the account of the Berry phase. It can be viewed as 
a probability of the condensate to change its momentum 
by 27r(po/2) due to an external potential. Indeed, by 
using a linear response relation Sp^^ ~ Xp^p [13]' where 
Xp is the static susceptibility of the electron gas, Pp^ is 
conventional (3D) conduction electron density and Vp is 
the disorder potential, one has (I^PpoP) = ^IXpoPl^oP' 
where A is the cross-section area of the wire. (Thus P3 
is proportional to the Born's scattering amplitude.) 

Since momentum pq is very high compared to both 
the Fermi momentum and the inverse interatomic dis- 
tance (recall that po — np^^A), we can expect that P3 is 
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rather smah. In order to estimate it, one needs to spec- 
ify Vp and Xp- At sufficiently high momenta the latter 
can be approximated by the susceptibility of a free elec- 
tron gas, Xp — Xp^oo — '^kp/{3TT-^Epp'^), where kp and 
Ep are Fermi momentum and energy (Ep = fi^fc|,/2m). 
For the disorder potential we assume that {V{r)V{r')) = 
V^^e"!''"'" I/", where the value of a is of the order of in- 
teratomic distance, a few A. The strength of the disor- 
der Vq can be estimated from the mean free path for- 
mula, vp/l = (Xp=o/4) / dn{l - cos6')(|Vp_p'p), where 
is angle between momenta p and p', whose magni- 
tudes are equal to kp, and Q is solid angle. Evaluating 
the Fourier transform of {V{r)V{r')) we have (|VpP} ~ 
Sttu^Vq / {1 + a^p'^)'^ , and, from the above mean free path 
formula we obtain that V^/a ~ {nkp)"^ /[2m'^nn{2akp)], 
where we assumed that {2akpY ^ 1- Then, we have 
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Note that po/kp is, up to a numerical factor, the number 
of transverse Fermi channels {N±) for the wire. 

It is instructive to estimate the QPS rate for a typical 
experimental system, such as 10 nm thick and 200 nm 
long MoGe nanowire. The value of conduction electron 
density p^^ — kp/Sn'^ can be obtained from the data on 
conductivity in the normal state. Using Drude formula 
a — e'^ p^^Tci/m, Td = l/vp, where the mean free path 
I is assumed to be of the order of interatomic distance, 
I ~ aX 14|, for a^^ ~ ..o^ f73i f,-^A h„ i /i-i 



we find kp '. 
p"^ ~ 3 X 10^** m-^ Then, for ^ = 8nrn _ 
a = AA we obtain that P3 ~ 10~^^. The phase ve- 
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locity Cs, which enters the first two factors in Eq. 
can be expressed in terms of the penetration depth A as 
Cs = cro/A(C)^/^ 0,11], where ro — 5nm is the cross 
section radius, c is the speed of light and the capacitance 
C is given by C ~ 27re[ln (L/rg)]"^, (e is the dielectric 
constant of surrounding medium). For a typical exper- 
imental situation e ~ 5 [l^ and so C ~ 9. For MoGe 
A ~ 0.72 [Til, which gives Cs ~ 7 x 10^ m/s and 
a ~ 2.4. Then we find that the attempt frequency, i.e., 
the first factor in Eq.® is loq ^ lO^^ - 10^^ m/s. The 
value of Sq decreases with the growth of the bias current. 
For / ^ 2fiA, which is of the order of critical current in a 
10 nm thick MoGe wire we obtain Sq ~ 3.2. Then, ac- 
cording to Eq.® Fqps does not exceed 10~^^s~^. Thus 
we conclude that the mechanism of condensate scattering 
at the disorder potential, e.g. Refs. 0,11, [ll[ is ineffective 
even in such disordered systems as MoGe nanowires. 

Note that the above estimate breaks down for the wires 
with the number of transverse channels N±_ of order 1. 
(For the above parameters N± ~ 100.) Also, in the pres- 
ence a weak link (a region with N± < 1 due to the 
inhomogeneity in wire's cross section, etc.) the main 
contribution into the integral in the QPS rate in Eq.® 
comes from the vicinity of point xnnk with /o(xiink) — 0. 



In that case one may set Ace ~ {xi, X2 — ^imk) in 
Eq. ^ and the QPS rate corresponds to that for a 
Josephson junction with an effective fugacity of the kinks 
- I / dxi exp [iiT J;^'^ dxp{x)]\. 

In the remaining part of the paper we investigate an- 
other source for the QPS production based on dissipa- 
tion of the plasmon mode due to the normal component 
of the electronic liquid. Such mechanism is motivated 
by the report of experimental observation of the QPS at 
high bias currents, i.e., when the value of the bias cur- 
rent is close to critical current @ . At such high currents 
superconductors can become gapless [16| , which leads to 
the appearance of appreciable normal component even at 
temperatures much lower than Tc. In order to estimate 
the effect of dissipation we recall that a superconducting 
wire can be described by an effective transmission line, 
which one may represent as a ladder of elementary build- 
ing blocks with inductors (or Josephson junctions) along 
one stringboard of the ladder and capacitors on the rungs. 
It is easy to show that in the continuous limit such model 
reproduces wave equation corresponding to the classical 
equation of motion for the Lagrangian in Eq. ®. The 
presence of the normal component can be modeled by 
adding resistors connected in parallel to the inductors. 
A straightforward analysis of the circuit shows that such 
resistors introduce dissipation according to the relation 
dE/dt = h^/{?>e^Rn) ^ dx{d^(t>/dTdxf, where i?„ is ef- 
fective resistance (per unit length) of the circuit and we 
have used the familiar relation between the voltage and 
the phase, AV^ — ftA(/)/2e. Such dissipative function cor- 
responds to an additional term in the action of the wire, 
which in the imaginary time representation can be cast 
in the form 
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where (J„ = l/RnA is the conductivity of the nor- 
mal component, an = e^p^^Ta\/m and is den- 
sity of normal electrons. Together with Cq this term 
yields familiar dispersion relation for the plasmon mode, 
a;2 - c2fc2 + iojclk^/ujd = 0, e.g. Ref. Q. 

The QPS rate in the presence of dissipation can be 
evaluated perturbatively in iSdiss- The first order term is 

rQPs = i^Im j d^T j dxidx2e-^^^^^'^Sai,,{dp,)/h, 

(12) 

where (/)c and S{(j)c) are given by Eqs. ® HJ, - here we 
neglect neglect the spatial variation of p, e.g., Eq. 
After a straightforward evaluation of integral in Eq. (fTTj) 
we obtain that 



^h^pj^'AT,, |Ax|3 
2m (Ax^ + c^Ar^ 



(13) 



Then, integration over At can again be performed with 
the use of Eqs. ® with an obvious replacement a 
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The first two factors in Eq. p4)) are identical (up to a 
numerical factor) to those in Eq. ([9|). The integral in 
the third factor can be evaluated by applying identity 
([6|). Note that unlike the previous case, e.g. Eqs. ([3] 
-0, the integrand function now contains an additional 
factor . Because the power of z is odd (thanks to the 
odd power of |a;| in 5diss in Eg. dill) ), the z-integration 
no longer produces biecspjl ± cos 6') (or its derivatives), 
which was the case in Eq. ^ and, as a result, led to the 
zero rate for CsP > //e. Evaluating the integral in the 
limit CsP I/e we obtain that the third factor (P3) in 
Eq. (HID 
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The rapid oscillations of the integrand function in Eg. p^ 
again lead to significant suppression of the rate. An 



estimate gives P3 
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where we have assumed 



again that I = 2jj,A, which is a typical value in exper- 
iments in Ref. [6[) and that ^ p^^ . This gives 
rqps ~ 10"'* . While this estimate is 4 orders of mag- 
nitude lower than the experimentally observed QPS rate 
0, it is 8 orders of magnitude greater than the previous 
one, e.g. Eqs.([9l [lOl). Moreover, in Eq.(fT4|) we did not 
account for the softening of the plasmon mode. Indeed, 
since Cg ^ (Ps)^^^ [2, sj, where ps is the density of super- 
conducting electrons, Cs is expected to drop to at the 
transition point (i. e., where pn = p and ps =0), which 
may lead to substantial increase in P3, etc. The discrep- 
ancy may also be related to other, more efficient sources 
of dissipation, such as shunting or contact resistances. 
A proper account of finite size of the nanowires may also 
lead to significant modifications of the QPS rates on both 
qualitative and quantitative levels [l9| . 

Finally we emphasize that the dissipation based mech- 
anism relies on the assumption that the fraction of the 
normal component is of order 1. Such an assumption 
obviously breaks down for a conventional (BCS) super- 
conductor, where at sufficiently low temperatures the 
density of the normal component is exponentially sup- 
pressed by the gap in the excitation spectrum. The 
normal component (and thus the relaxation rate of the 
plasmon mode) can be significantly enhanced by pair- 
breaking mechanisms, such as high bias current (as in 
Ref. [6|). Note that the plasmon resonances in super- 
conducting nanowires can be studied experimentally by 
measuring the reflectivity coefficients of the wires, e.g., 
Ref. (2^. Such measurements could therefore test the 
validity of arguments presented in this Letter. 



In summary we have shown that a proper account of 
the momentum conservation during the QPS formation 
leads to a very strong suppression of the QPS rate, e.g. 
by factor P3 in Eqs.([9l ITOl) (for the potential scattering 
mechanism) and P3 in Eqs. ([131 [T5t (for the dissipation 
based mechanism). The latter mechanism, under spe- 
cial circumstances, may lead to an observable QPS rate 
and may be the origin of the enhancement of phase slip 
transitions observed in Ref. 6] . 
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